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f^*) , Abstract 

■ The eigenvalue correlations of random matrices from the Jacobi Unitary Ensemble have 

£N) | a known asymptotic behavior as their size tends to infinity. In the bulk of the spectrum 

the behavior is described in terms of the sine kernel, and at the edge in terms of the Bessel 
^ (— i kernel. We will prove that this behavior persists for the Modified Jacobi Unitary Ensemble. 

O This generalization of the Jacobi Unitary Ensemble is associated with the modified Jacobi 

weight w(x) = (1 — x) a (l + x) /3 h(x) where the extra factor h is assumed to be real analytic 
and strictly positive on [—1,1]. We use the connection with the orthogonal polynomials with 
respect to the modified Jacobi weight, and recent results on strong asymptotics derived by 
K.T-R McLaughlin, W. Van Assche and the authors. 



X 1 1 Introduction 



In the early sixties, Dyson predicted that the local correlations between the eigenvalues of 
ensembles of random matrices, when their size tends to infinity, have universal behavior in the 
bulk of the spectrum. He expected that this universal behavior depends only on the type of 
the ensemble: orthogonal, unitary or symplectic. This constitutes the famous conjecture of 
universality in the theory of random matrices. For the classical ensembles (Hermite, Laguerre 
and Jacobi), this conjecture has been proven, see for example [[0| 18, [l9|, 23]. For the unitary 



ensembles much more is known due to the connection with orthogonal polynomials, and the 
universality conjecture in the bulk of the spectrum is proved for a wide class of unitary ensembles, 
see |, |, lH. 



At the edge of the spectrum this universal behavior breaks down. For Hermite ensembles, it 
is known that the local correlations (at the soft edge) can be expressed in terms of Airy functions 
[ 2], [lC|, |27|] , and for Jacobi and Laguerre ensembles (at the hard edge) in terms of Bessel functions 
10, 17^2^, 28 1 . For example, for the Jacobi Unitary Ensemble 
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where w(x) = (1 — x) a (l + x)@ is the Jacobi weight, the eigenvalue correlations near 1 are 
expressed in terms of the Bessel kernel 

Jf (y, v) = Jai^VvJ'atVv) ~ Ja(Vv)VuJ'a{Vu) /j g N 

2(u — v) 

as n — > oo. J a is the usual Bessel function of the first kind and order a. The order agrees with 
the exponent of 1 — x in the Jacobi weight. 

Nagao and Wadati [^0], §6] expect that a universality result persists for more general Jacobi- 
like ensembles, in the sense that the local form of the weight function near 1 determines the 
eigenvalue correlation near 1. It is the aim of this paper to prove this universal behavior for 
a generalization of the Jacobi Unitary Ensemble, which we call the Modified Jacobi Unitary 
Ensemble (MJUE). The MJUE is given by ( [1,11 ) with modified Jacobi weight 

w(x) = (1 - x) a (l + xfh(x), for x £ (-1,1), (1.3) 

where a,f3> —1 and the extra factor h is real analytic and strictly positive on [—1,1]. The 
Modified Jacobi Ensemble is a probability measure on the space ofnxn Hermitian matrices 
with all eigenvalues in (—1,1). The MJUE gives rise to a probability density function of the n 
eigenvalues x±,X2, ■ ■ ■ ,x n given by 



1 - 

P (n \x 1 ,x 2 ,... ,x n ) = -^-[[wix^lllxi ~Xj\ 2 , (1.4) 

n 3=1 i<3 

with xi, X2, ... , x n £ (—1, 1) and Z n a normalizing constant. 



T^n,m{xi , • • • 3?m) — / -. , / • • • / P (^1 j • • • j %m > 2-m+l j • • • ; %'n)dXm+l • • • dx r 



Dyson || showed, see also 16], that we can express the correlation functions 7^ n ^ r 

n\ 

(n — m) 

in terms of orthogonal polynomials. Denote the nth degree orthonormal polynomial with re- 
spect to the modified Jacobi weight w by p n (z) = p n (z;w) = r ) n z n + "f n > 0. Then 
R-n,m{xi, ■ ■ ■ x m ) = det (K n (xi, a;j))i<ij< m > where 

n-l 

K n (x,y) = y / w(x)y / w{y)^2p j {x)p j (y). (1.5) 

j=o 

By the Christoffel-Darboux formula, we have 

K n (x,y) = \/w(x)yJw{y) , (1.6) 

In x-y 

which shows that asymptotic properties of K n are intimately related with asymptotics of the 
orthogonal polynomials p n as n — > oo. 

In a previous paper with K.T-R McLaughlin and W. Van Assche [12|, we studied the asymp- 
totics of the polynomials that are orthogonal with respect to the modified Jacobi weight. We 
used the Riemann-Hilbert formulation for orthogonal polynomials of Fokas, Its, and Kitaev |9 
and the steepest descent method for Riemann-Hilbert problems of Deift and Zhou [0]. In []12 
we concentrated on the asymptotics of the polynomials away from the interval [—1,1], but the 
Riemann-Hilbert method gives uniform asymptotics in all regions in the complex plane. Here 
we are interested in the behavior on [—1,1], and in particular near the endpoints ±1. The 
Riemann-Hilbert method was applied before to orthogonal polynomials by Deift and co-authors 
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|3|, H|, ||, |ll]]. They studied orthogonal polynomials on the real line with varying weights, and 
used the asymptotics to prove the universality in the bulk of the spectrum for the associated 
unitary ensembles. We apply the same method to prove the universality at the edge of the 
spectrum for the MJUE. Our main result is the following. 



Theorem 1.1 Let w be the modified Jacobi weight ( |l.q ) and let K n be the kernel fll.q ) associated 
with w. Then the following holds. 

(a) For x G (—1, 1), we have osn-> oo, 

-K n (x, x) = 7 =L= + o(-). (1-7) 

The error term is uniform for x in compact subsets of (— 1, 1). 

(b) Let £(x) = =g . Then for x G (—1, 1) and u, v G E ; we /iaue as n —* oo, 

1 / it t> \ sin7r(ti — u) „/l\ , 
-tf„ x+— — ,a; + — — = j-i -L + ol-). (1.8) 



n£(x) \ n£(x)' n£(x) J ir(u — v) \n 

The error term is uniform for x in compact subsets of ( — 1,1) and for u,v in compact 
subsets ofR. 

(c) For u, v G (0, oo), we have as n — > oo, 

i £ »( 1 -i' 1 -i)= , - ( "-"» + °(^)' (L9) 

where J a is i/ie Bessel kernel given by (|1.2|). T/ie error term is uniform for u, v in bounded 
subsets of (0, oo). 

Note that the error term in ( |1.9| ) holds uniformly for u,v in bounded subsets of (0, oo), not just 
in compact subsets. By symmetry, there is a corresponding universality result near —1. 

The eigenvalue density is the 1-point correlation function TZ n ^{x) = K n (x,x), see for 
example |l6| . Therefore, part (a) of the theorem yields the asymptotic eigenvalue density 



T^-n,i{ x ) ~ n£(x) as n — > oo. This result is in agreement with jl3|, |1S|| . The scaling in (l.S) 
has the effect that x is the new origin and that the asymptotic eigenvalue density at x is 1. 
At the endpoints, (1.7) breaks down, and the eigenvalue density is 0(n 2 ) as n — > oo, near the 



endpoints, see for example |13| . This explains the scaling in (|L 

Part (b) of the theorem states the universality (independent of the choice of a, (3, h and x) 
for K n in the bulk of the spectrum. It extends the result of Nagao and Wadati Jl9|, (4.19)] for the 
case that h = 1. At the edge 1 of the spectrum we have a universality class for K n (independent 
of the choice of (3 and h) which is only affected by the local form of the modified Jacobi weight 
near 1, see part (c). 

Using Theorem [O] we can answer local statistical quantities concerning the eigenvalues. 
Here we follow || |j. The probability P n (a,b) that there are no eigenvalues in the interval 
(a, 6) C (—1, 1) is given by 

P n (a,b) = det(I - K n ), 

where K n is the trace class operator with integral kernel K n (x,y) acting on L 2 (a, 6), and where 
det(I — K n ) is the Fredholm determinant. For a fixed interval (a, 6) we have that P n (a, b) — > 0, 
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as n — > oo. So, to understand the asymptotic behavior of P n at the edge of the spectrum, we will 
look at intervals near the edges which shrink with n, and we are led to consider the asymptotic 
behavior of P n (l — 7^7, l) as n — > 00, where s > 0. We have the following universality for P n 
at the edge 1 of the spectrum, depending on the parameter a but independent of the choice of 
(3 and h. 

Corollary 1.2 For s > 0, we have 

lim P n (l - -£-, l) = det(J " 3a,s), (1.10) 

where is i/te integral operator with kernel § a (u,v) acting on L 2 (0,s), and det(J — S a ,s) 
the Fredholm determinant. 

As mentioned before, our main tool in proving Theorem 1.1 is the asymptotic analysis of the 
Riemann-Hilbert problem for the orthogonal polynomials with respect to the modified Jacobi 
weight, as developed in jl^]. We give an overview of this work in section 2. This approach is 
able to give strong and uniform asymptotics for the orthogonal polynomials in every region in 
the complex plane, which we also review in section 2. The proofs of Theorem 1.1 and Corollary 
1.2 are given in section 3. 



2 Riemann-Hilbert problem for orthogonal polynomials 



In this section we will recall the Riemann-Hilbert problem (RH problem) from E^] for the 
orthogonal polynomials for the modified Jacobi weight, and the steps in the steepest descent 
method that are used to obtain the asymptotics on the interval. 

The analysis starts from a characterization of the monic orthogonal polynomials 7r n for the 
modified Jacobi weight w given by (|b^) as a solution of a RH problem for a 2 x 2 matrix valued 
function Y{z) = Y(z;n,w). This characterization of orthogonal polynomials is due to Fokas, 
Its, and Kitaev ||. The conditions ( [2.3; ) and (2A) are needed to control the behavior near the 
endpoints, see [[[^] for discussion. 

(a) Y(z) is analytic for z G C \ [—1,1]. 

(b) Y possesses continuous boundary values for x G (—1,1) denoted by Y + (x) and Y^(x), 
where Y + {x) and Y_(x) denote the limiting values of Y(z) as z approaches x from above 
and below, respectively, and 



Y + (x) =K.(x) 



w(x) 
1 



for x £ (—1, 1). 



(2.1) 



(c) Y{z) has the following asymptotic behavior at infinity: 



Y(z) =[1 + 



z n 
z~ n I ' 



as z — > 00. 



(2.2) 



(d) Y(z) has the following behavior near 2 = 1: 

1 |z-l| Q 



Y{z) = I 



O 
O 



1 |z-l |Q - 

1 log 1 2 — 1 
1 log 1 2 — 1 
1 1 



o 



1 1 



if a < 0, 
if a = 0, 
if a > 0, 



(2.3) 



as z -► 1, z G C\ [-1,1] 
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(e) Y{z) has the following behavior near z = 



-1: 



Y(z) 



O 

o 



1 \z + lf 
1 \z + lf 
1 log \z + 1 
1 log I z + 1 

1 1 
1 1 



O 



if p < 0, 
if /? = 0, 
if /? > 0, 



(2.4) 



as z 



■MeC\[-i,i]. 

The unique solution of this RH Problem is given by 



(2.5) 



where 7r n is the monic polynomial of degree n orthogonal with respect to the weight w and with 
7 n the leading coefficient of the orthonormal polynomial p n . 

We apply a number of transformations Y*-^T>-^St-^Rto the original RH problem in 
order to arrive at a RH problem for R, which is normalized at infinity, and whose jump matrices 
are uniformly close to the identity matrix. Then, R is uniformly close to the identity matrix, 
and by tracing back the steps we deduce the asymptotic behavior of Y. 

In the first transformation we turn the original RH problem into an equivalent RH problem 
for T, which is normalized at infinity and with a jump matrix whose diagonal elements are 
oscillatory. Let 



0-3 



1 

-1 



be the Pauli matrix, and let T be given by 
T(z) = 2 na3 Y{z)^{z)- na \ 



(2.6) 



where ip(z) 



z + (z* 



1)V2 for zE C\[-l,l], so that (f is the conformal map from C \ [— 1, 1] 



onto the exterior of the unit circle. Then T satisfies a RH problem which is normalized at infinity 
(i.e., T(z) — > I as z — > oo), and 



T, (x) = T(x) 



<p+(x) 2n w(x) 

ip_{x)- 2n 1 ' 



for x £ (—1,1). 



-1 



Figure 1: The lens S 
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The jump matrix for T factors into a product of three matrices. Using this factorization, 
we transform the RH problem for T into an equivalent RH problem for S, with jumps on a lens 
shaped contour S, as in Figure 0. S is defined in terms of T by 



T(z), for z outside the lens, 



T ( z ) (- w (z)-^cp(z)- 2n l) ' f ° r Z in the Upp6r Part ° f the lenS ' ( 2 - 7 ) 

T(z) ( , , i , , 9 „ „ | , for z in the lower part of the lens. 
v ' \w(z) (p{z) zn 1 1 



S(z) 



For the third transformation, we need to construct parametrices in the outside region and 
near the endpoints ±1. Constructing the parametrix in the outside region, we need the Szego 
function D(z) = D(z; w) associated with the weight w given by 

. . (iz 2 -!) 1 / 2 f 1 logw(x) dx \ r r n 

D(z) = exp 1 - > / t—4 , for z£C\[-l,l]. (2.8) 



2tt J-i^/T^z-x j' 

The Szego function is a non-zero analytic function on C \ [—1, 1] such that D + {x)D^.{x) = w(x) 
for x £ (—1,1). The parametrix N in the outside region is given by 

(a(z)+a(z)~ L a(z)—a(z)~ 1 
,(,)- 2 aW -i (,)4)-^ i)UrT! (L,,)! 
-2i 2 

where = lim^oo D(z) and a(z) = (z - l) l /*(z + l)" 1 / 4 . 

Next, we define a parametrix P in Us, which is the disk with radius 5 and center 1, where 
5 > is sufficiently small, by 

P{z) = E n (z)y(n 2 f(z))W(z)-^<p(z)- n ^, (2.10) 

where f(z) and W(z) are scalar functions given by 

/(;.) = ^ (log <^)) 2 , (2.11) 

and 

W{z) = (z- l) a l 2 {z + \fl 2 h x l 2 (z). (2.12) 

In ( |2,10p ^(C) is a 2 x 2 matrix valued function defined for £ 6 C \ where is the 
contour shown in Figure ||. For our purpose here, it suffices to know the expression of ^(C) 
for 27r/3 < arg£ < tt. This is given in terms of the Hankel functions Ha^ and • For 
27r/3 < argC < tt, we have 

_/ ifli 1 »(2(-C) 1 /=) \hT\2(-Q^) \ 1 . 

* (c) = (rf»)'( 2 (-o ,/2 ) < 1/2 Ww-o^J (2 ' 13) 



The factor E n (z) in (l2.1C| ) is given by 



£ n (z) = iV(z)iy(z)-3_L ^. M /(zr /4 (27rnr3 /2 (2 M) 

E n (z) is analytic in a full neighborhood of Us- 
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2ji 

3 



Figure 2: The contour 



There is a similar definition for the parametrix P in a 5 neighborhood U$ of —1, see [12] for 
details. 

We then have all the ingredients for the third transformation. We define 
R(z) = S(z)N- 1 (z), for z€C\ (U s (JU 5 UE), (2.15) 
R(z) = S(z)p- 1 (z), iovzeU s \E, (2.16) 
R{z) = S(z)p- 1 (z), iovzeU s \E. (2.17) 
Then R satisfies a RH problem with jumps on the system of contours Sjj shown in Figure 0. 





Figure 3: System of contours 

Note that R depends on n, but the contour Y,r does not depend on n. The jump matrices 
for the RH problem for R turn out to be uniformly close to the identity matrix with error term 
O(^). Then it follows that 

R(z)=I + o(^\, asn^oo, (2.18) 
uniformly for z G C \ T,r, and also, see || Section 8.1] in particular formulas (8.19) and (8.20), 
-^R(z) = 0(1), asn^cxD, (2.19) 

uniformly for z £ C \ Ejj. In the following we will also use that 

detR(z) = l. (2.20) 



Remark 2.1 The asymptotic result ( 2.18; ) for R may be refined and it is possible to obtain a 
full asymptotic expansion 

R(z) ~ I + — — -r— , as n — > oo. 
^-^ n K 
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The matrix coefficients Rk{z) are explicitly computable in function of a, (3, and the analytic 
factor h in the modified Jacobi weight. The first two coefficients are computed in Jl~2| ]. 



Remark 2.2 In [|T^] we also derived asymptotic expansions for the polynomials p n and 7r n 
uniformly valid in compact subsets of C\ [— 1, 1], as well as asymptotic expansions for the leading 
coefficient 7 n and for the coefficients in the recurrence relation satisfied by the orthonormal 
polynomials. For example, we have by [12, Theorem 1.4] uniformly for z E C \ [— 1, 1], 



2 n ir n {z) 
ip(z) 



<p(z 



,1/2 



D(z) ^2>2 _ !)i/4 



1 + 



as n 



oo. 



n 



The 0(£) term can be developed into a complete asymptotic expansion in powers of 

The Riemann-Hilbert method also leads to strong asymptotics on the interval (—1, 1) and 
near the endpoints dbl. While these results are closely related to the asymptotics of K n as given 
in Theorem 1.1, we do not actually rely on them in the proof of Theorem 1.1. Therefore we 
state here the asymptotics of the orthogonal polynomials without proof. The proof is similar to 
the proof of Theorem 1.1, and in fact somewhat simpler. See also pj, |!|. 

As before, we let 5 > be the radius of the disks Us and Ug. Then we have for x £ 
(-1 + 6,1-6), 



7T„(x) 



V2D Q 



2™ v /^x)(l - x 2 ) 1 / 4 



Rn(x) cos ( (n + 1/2) arccos x + ip(x) 



7T 



D 2 



R\2(x) cos ( (n — 1/2) arccos x + ip(x) 



(2.21) 



where Rn(x) = 1 + O(^) and Rn(x) = O(^) as n — > oo, uniformly for x G (—1 + 6, 1 — 6). 
The O(^) terms have a complete asymptotic expansion, see |l^]. The function ip(x) in ( 2.21|) is 
given by 



2p(x) 



vl — x 2 l' 1 \ogw(t) dt 



2ir 



£ 



a(arccos x — it) + (3 arccos x I + 



Vl-x 2 f 1 log h(t) dt 



2tt 



V\—pt-x' 



(2.22) 



where the integral is a Cauchy principal value integral. We remark that, under various assump- 
tions, asymptotic results on the interval of orthogonality have been established by many authors, 
see for example (§ § 0, [H| g§ |§. 
For x € (1 — 6, 1), the following is valid 



7T n (x) 



n arccos x 



,V2 



2«v^) (1-X 2 )V4 

[i?ii(x) (cos (i(x)J a (n arccos x) + sin d(x)J' a (n arccos x)) 

— — 2~Ri2(x) (cos(2(x)J a (naicccosx) + sin C t 2{x)J l ct (n arccos a;)) 

where R\i{x) = l + 0(~) and #12 (a) = 0(=r) 

as n — ► oo, uniformly for x G (1 — (5, 1), and where 
1 



(2.23) 



07T 



arccos x + ^>(x) + r> . 



Ci(x) = - arccos x + ^(x) + — , ( 2 (x) 
with ijj given by ( |2.22j ). There is an analogous expression for x in the interval (—1,-1 + 5). 
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3 Proof of Theorem 1.1 and Corollary 1.2 



In this section we prove Theorem 1 . 1 and Corollary |1.2| . We follow the work of Deift et al. (3|, Q . 
Recall that 



In x-y 

= V w (x)V w (y) r yn-i • 



x-y 

where ir n is the monic orthogonal polynomial of degree n with respect to the modified Jacobi 
weight w. As in (3|, || we replace the polynomials 7r n _i and 7r n by the appropriate entries of Y, 
see (12.51), to obtain 



z, r s 1 r~r^ r ^ Y ll {x)Y 2l {y) -Y 2l (x)Y ll {y) 

Kn{x,y) = - — V u>{x)V w (y) • 

2ni x — y 

Thus, K n can be expressed in terms of the first column of Y. The asymptotic behavior of Y 
follows from the transformations Y i— > T i— ► 5 i— ► i? described in section 2, and the behavior 
pl8|)-(|2T20D of i?. 



3.1 Proof of Theorem 1.1 (a) 

We will first express Y\\ and Y21 in terms of R. In the following, 5 > will be a small but fixed 

number. This number is the radius of the disks Us and Ug used in the local RH analysis around 
±1. 

Lemma 3.1 We have for x G (— 1 + S, 1 — S), 

/ Yn(x)\ 1 ncr ( e inwccosx \ 

v (^\ ) = / / \ ^ 3 -L+(x) I i narccosa; I 1 (3-2) 
with L(z) given by 

( a(z)+a(z)- 1 a(z)-a(z)- 1 \ 

^ 2 ^ 1 ^ 2 ^ ! e ^ )<T3 > ( 3 - 3 ) 
-2i 2 / 

u;aere 

wia;) = - afarccos x — tt) + p arccosa; H + , . 3.4) 

V 7 2V v ' ) 2vr J^^/rrpt-x K ' 

The matrices L + (x) and -^L + (x) are uniformly bounded for x G (—1 + 5, 1 — 5) as n — > 00, and 

detL + (x) = l, /or x e (-1 + 0, 1 -0). (3.5) 



Remark 3.2 The singular integr al fl3.4j ) is being understood in the sense of the principal value. 
It may be shown that D + (x) = w{x) ex.p(—iip(x)) so that —ip(x) is the argument of the Szego 
function on the interval. 

Proof. We use the series of transformations Y T S 1— > R and we unfold them for z in the 
upper part of the lens but outside the disks Us and Us, and then take the limit to the interval 
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(— 1 + 5,1 — 5). Thus, let z be in the upper part of the lens but outside the disks Us and Us- We 
then have by Q, Q and (pl^) 

y W = 2— jk, W ,) ( ( /_ ( $ z) _„ ^ J ■ (3.6) 



Inserting the expression ( |2.9| ) for JV into ( |3,6| ), we obtain for the first column of Y 

- 2— »J2(«)I5g _ * _ ,* J/ I • (3-7) 



We now take the limit z — ► x S (—1 + <5, 1 — <5). Since -D+(x) = \J w(x)e~ v ^^ and ^+0*0 



exp(i arccos x), (|3.2[) now follows from (3.7), 

Note that i?(x) and -^R(x) are uniformly bounded for a; G (—1 + 5, 1 — <5) as n — > oo. Let J7 
be a neighborhood of [—1, 1] such that log/i is defined and analytic in U, and let 7 be a closed 
contour in U\ [—1,1] encircling the interval [—1,1] once in the positive direction. Via contour 
deformation, we may write tp in the form 



, . 1 / / . _ \ vi -x 2 r iog/i(c) d( 

ip(x) = — (^a(arccos x — ir) + p arccos xj H — : — J> 



( C 2_ l) l/2 C _ X - 



Thus ip has an analytic extension to a neighborhood of (—1,1). This implies that ip and its 
derivative are bounded on [—1 + 5,1 — 5). From the explicit form (|3.3| ) of we then find that 



L + (x) and -^L + (x) are uniformly bounded for x S (—1 + 5, 1 — 5) as n — ► 00. Since det-R(-z) = 1 
it is easy to see from (|3.3[) that detL + (x) = 1. □ 



Proof of Theorem 1.1 (a) Letting y — ► x in ( |3.1| ) we get 

k-„(x,x) = ^(x)(y 11 (x)y 2 ' 1 (x)-y 21 (x)y 1 , 1 (x)) 



— det rvM^^i^) vMxjy/^x) 

2vri Vv^50 y 2l(a:) TMxT^Vx) 
1 ( 2*y/Mx)Yu{x) ±(2 n ^(x)Y 11 {x)) 

— det / \ . (3.c 

2""v^(xTy 21 (x) £ (2-™ v ^(xTy 21 (x) x 



2vri 



By (ET^) the matrix in (ET3) is equal to 



'in arccos 



„m arccos x in „in arccos x \ / ; \ /~ 

] + f-L + (x)V° 6 

g— inarccosz in Q~ ln arccos x j ' I dx j \0 e — marccosa; 



n in arccos x 



L , (xW narccosx<J3 Vi-x' 2 1 I —L,(x) 

+ { ' I ^ m J ~ \^dx J \0 e~ l " am '" s ■'' 

Since L+(x) and e m arccos xo " 3 have determinant one, it follows that 



K n (x,x) = —— det 
2m 



1 — \ / <7 \ /n w 

^ vT^ + e - 4 n arccos ^3 L -l (x) ( ^-L+{x W 



in arccos x 



q g—«n arccos x 

(3.9) 
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The entries of L + (x) and -jj^L + (x) are uniformly bounded by Lemma 3A. Since detL + (x) = 1, 
also the entries of L^}(x) are uniformly bounded. Thus, we have uniformly for x £ (—1 + 5,1 — 5), 



K n (x,x) 



1 



n 



7T\/l — x' 



VT- 
+ 0(1). 



2in 



2™ U 77T=5 + 0(1) 2-ai\y/l^ 



+ 0(1) 



This proves part (a) of Theorem 1.1. 
3.2 Proof of Theorem 1.1 (b) 



□ 



and v x>n to denote x + n ^ x ^ and x + n ^ x ^ , respectively. We write 



"K n (u x n , V x n) 



Proof of Theorem 1.1 (b) Let x £ (—1, 1) and u, v G R. For the sake of brevity, we use u x ^ n 
X} n to dei 

K n (u,v) - — 
n^(x ) 

then have 

Y\i (",-.») Y\i( )',-.„) ^ 

1 \ (3.10) 



From (3.1) we then have 



K n (u,v) = -7—7 rJw(u xn )Jw(v xn ) det ■ , . , 

2m(u-v)y 'V \Y 2 i(u x ^) Y 2 i(v x>r 



1 



det 



2" v / w(u Xin )Y n(u :Ein ) 2 n v / w(z; a;iTO )y ii(?; a;ira ) 

2m(u - V) V 2 ~ n V / H^^ y 2l(^,n) ^~ n y/'w(Vx,n)Y2l(v X) n 

Next, we use the expression (|3.2| ) for Y"u and Y21 in ( |3.10| ) to obtain 



K n (u,v) 



2m(u — v) 
+ L + (v x>n ) 



1 

2ni(u — v) 



det 



L + (u x ^ n ) 



„m arccos u x 





g— in arccos u x , n Q 



e 
e 



in arccos v x n 



■in arccos v x n 



det 



L+(v x>n ) 



+ [L + (u x n ) - L + (v x<n )) 



m arccos u x , n 



-111 arccos u x 



in arccos u x n 



in arccos Vx.n 



-in arccos v x n 





g— in arccos u x , n q 



(3.11) 



Since det L + (v x ^ n ) = 1, we then get 



K n (u,v) 



1 



2m(u — v) 



det 



m arccos 



-in arccos u x 



+ L, 1 (v X)n ) [L + (u x>n ) - L + (v xn )] 



m arccos v x n 



-in arccos v x 



.-,111 arccos u r 



e -in arccos u x ,n Q 

Since -^L + (x) is uniformly bounded asn-> 00, we have by the mean value theorem 







(3.12) 



L + (u Xtn ) - L + (v x<n ) = O (u, 



x,n u x,nl 



o 



u — v 



11 



uniformly for x £ (—1 + 5,1 — 5) and for u, v in compact subsets of M. Since L + (x) is uniformly 
bounded, and detL + (x) = 1, we have that L~7 (v x ,n) is uniformly bounded as well, so that 



x,n l\ 



m arccos u x ,n 



_/o(V) ° N 

e ~m arccos u XtTl q J I q ^ u -v^ q 



(3.13) 
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Thus by and 



^ Y ( e in a,rccosu Xt n _|_ Q f u—v ^ g 

K n (u,v) = —-z-T-, rdet , 71 . 

! 1 -"> " vm 9%,n _|_ Q ntz2L\ g" 



in arccos i.> ;T 



2vu 


;«-«) 




i 


27ri 


[u-v) 




1 



«n arccos v x ,n 



c vn arccos u x n c m arccos v x n \ / i 



27TZ(W 

sin(n(arccos i^.n - arccos v X)U )) + Q / 1 ' , . ( _ 



7r(n — f ) \n 
uniformly for x £ (—1 + 5, 1 — 5) and u, v in compact subsets of M. Since 

n(arccos u Xjn — arccos v X)n ) = —ir(u — v)(l + 0(l/n)), as n — > oo, 

also uniformly for x £ (—1 + 5,1 — 5) and u, t> in compact subsets of M, part (b) of Theorem 1.1 
follows. □ 

3.3 Proof of Theorem 1.1 (c) 



To prove part (c) of Theorem 1.1, we start with a result similar to Lemma 3A_ 
Lemma 3.3 For x € (1 — 5, 1), we have 



Yii(x)\ _ I 2vrn o _ n(T3 n/f ( ^( J Q (narccosx 



v i2l(ic)y y w(x) \ § arccos xJ' a {n arccos x 

mf/i M(z) given by 

1 f 1 ~A /Y*W 4 



M(z) = i?(^)iV^)W(z) CT3 -^ I _. ^1 /(z)^ 4 , (3.16) 

where R is the result of the transformations Y T ^ S i— > i? o/ f/ie problem, the matrix 
valued function N is given by (2.9), and t/ie scalar functions f and W are given by ( |2.11 ) and 



( 2.12 ), respectively 

lit 



M is analytic in Us with M(z) and 4-M{z) uniformly bounded for z £ Us as n — > cxd. 



Furthermore, we have 
detAffjs) = 1. 



Proof. As in the proof of Lemma 3.1 we unravel the transformations Y i— ► T ^ S *—> R, but 
now for z in the upper part of the lens and inside the disk Us- We then have by (E^6j), (f2~7|), 
(|2l0D and ( gig ) 



K(z) = 2-~»Riz)E n {z)9 (n 2 f(z)) W(z)— fj . (3.17) 



Since Imz > 0, we have by ( 2.12j) that VF(#) = w 1 ^ 2 {z)e ma ^ 2 . Inserting this into ( 3.17 ) we get 
for the first column of Y 



(£(*)) =^ 1/2 (^)2- nCT3 ^)^(.)^(n 2 /(^)) ( 



g— iwa/2 

e 



ina/2 J ' ( 3 - 18 ) 
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Since z is in the upper part of the lens and inside the disk Us, we have 2tt/3 < aign 2 f(z) < n, 
see [E|, section 6], and we thus use ( p,13 ) to evaluate ^(n 2 f(z)). From formulas 9.1.3 and 9.1.4 
of ll I we then have 



-iira/2 

* (n 2 f{z)) ( e ^ a/2 



(3.19) 



J a (2n(-/(z)) 1 /2) \ 
^nfy\z)J' a {2n{-f{z)f/ 2 )) ' 

By (pHl ) and ( jOBD we have R(z)E n (z) = M(z)(2vrn) CT3 / 2 . Inserting this and ( |3~!9D into ( glf ) 

we get 

J Q (2n(-/(^))V2) \ 



(3.20) 



We now take the limit z — > x G (1 — (5,1). By (2.11), and since 92+(x) = exp(i arccos x) we 
have /^ 2 (x) = | arccosx, so that 2(— /^(x)) 1 / 2 = arccosx. Inserting this into fl3.20|) , we obtain 



M is analytic in [7,5 since both £7 n , see [12, Proposition 6.5], and i? are analytic in So, 
we may write M(x) instead of M + (x) in ( |3.15| ). 

By ( 2.18| ) and (|2.19j ) we have that R{z) and -^R(z) are uniformly bounded for z £ [/j as 
n — > oo. Since 



0-3/4 



is analytic for z G C/5 and does not depend on n, we have from flOp that also M(z) and ^M(z) 
are uniformly bounded for z G C/^ as n — > 00. 

The fact that detM(z) = 1 follows easily from (pl| , (|2T20| ) , and flOf ). □ 



Remark 3.4 We check that ( gig ) yields ( y 11 y J = 0(1) /" 1, which is in agreement 

with the fact that Y"n and I21 are polynomials, see (2.5). Since arccosx = O ((1 — x) 1 / 2 ) as 
x /" 1, we have by formula 9.1.10 of O 



J a (n arccos x) 
arccos x J' a (n arccos x^ 



C>((1 

o((i 



x) a/2> 
:rW 2< 



as 1 / 1. 



(3.21) 



Since M is analytic near 1, see Lemma |3.3| , we have M(x) = O(l) as x / 1. Inserting (|3.21j) 
into ( [3.15 ) and noting that w(x) _1 = 0((1 — x)~ a ) as 1; / 1, we then have indeed that Yn(x) 
and ^2i(x) remain bounded as x /* 1. 

We also need the asymptotic behavior of J a (u n ) and J' a {u n ) as n -> 00, where we put 
u n = 1 — 2^ and u n = n arccos u n . These will be contained in the next lemma. 



Lemma 3.5 Let u G (0, oo), u n = 1 



' 3 
1i 2 



71- 



2^-, and tt n = n arccos u n . We then have as n — > 00, 

(3.22) 



Ja(u n ) = Ja(Vu) + O 



j'a(Un) = JaiVu) + O 



-+1 



71- 



77ie error terms hold uniformly for u in bounded subsets of (0, 00) 



(3.23) 
(3.24) 
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Proof. Since z 1//2 arccos(l — z) is analytic in a neighborhood of with expansion 



-1/2 



arccos 



1 - z) = V2 + 0(z), as z ^ 0, 



we easily get fl3.22|) . 

By formula 9.1.10 of |l]] we know that J a (z) = z a G(z) with G an entire function. From 
( 3.22 ) and Taylor's formula we then get uniformly for u in bounded subsets of (0, oo), 



J a (n arccos u Tl 



u a 1 + O 



ii- 



G(Vu) + O 



' 1 



11, 



-+1 



which is ( |3.23|) . The proof of fl3.24|) follows in a similar fashion, using 

J' a (z) = az a ~ x G{z) + z a G'(z), 
and Taylor's formula again. 

Now we are ready for the proof of our main result. 
Proof of Theorem 1.1 (c) Let u,v £ (0, oo) and define 



□ 



1 



2n 2 ' 



2n 2 ' 



We put 



u n = n arccos u n , v n = n arccos v r . 



(3.25) 



D n (u,v) = -^K n (u n ,v n ) 



From ( |3.l|) we then have 



n i \ 1 r~t — V r~i — \a + f Y n(u n ) Yxi(v n ) 

D n {u,v) = — — ry/w{u n )y/w{v n )det [ , ' , ' 

2m(u-v) v /v \Y 2 i{u n ) Y 2 i{v n ) 



2" v / ^K)y 11 (u w ) 2 n y/ M^)Y u{v n ) 

27ri(u - v) ™" V2-"v / ^K) y 2iK) 2-™vVK)y 21 K) 



■det 



(3.26) 



(3.27) 



Next, we replace the two columns in the determinant in ( |3.27 ) by the expression ( 0.15 ) we found 
in Lemma |3.3|. It follows that 



D n (u,v) 



2irn 



2iri(u — v) 



det 



M(u n ) 



J a (u n ) 



+ M(v n ) 



J a (v n ) 
2n 



hv n J'Jy n ) i 



(3.28) 



We rewrite the matrix appearing in the determinant in ( |3.28| ) as 



M(v n ) 



+ M(v n )- l [M(u n ) - M{v n )} 



Ja(u n ) N 

^u n J' Q (u n ) 



Now we use detM(u n ) = 1 and the fact that M(z) is uniformly bounded for z £ Us, see Lemma 
3.3 , to conclude that the entries of M(v n )~ l are uniformly bounded. By Lemma 3^, we also 



have that -^M(z) is uniformly bounded so that M(u n ) — M(v n ) = O (^f)- From Lemma 3.5 
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it follows that J a (u n ) = 0{u a / 2 ) and u n J' a (u n ) = 0(u a l 2 ) uniformly for u in bounded subsets 
of (0, oo) as n — > oo. Hence we have, uniformly for u, v in bounded subsets of (0, oo), 

M{v n )~ l [M{u n ) - M(v n )] V / = ) n J . (3.29) 

It now follows that (we use detM(u n ) = 1) 

1 / J a (u n ) + ofeu§) J a (v n 

D n {u,v) = - -det { i\ h ( 3 - 3 °) 

2 ( U ~V) \u n J>(u n )+0(^Ul) V n J' a (Vn) 



Since J a (v n ) = 0(v a / 2 ) and v n J' a (v n ) = 0(v a ^ 2 ) as n —* oo, we then get uniformly for u,v in 
bounded subsets of (0,oo), 

1 / Ja(u n ) Ja(v n ) \ ( U ^V^\ 

D n (u,v) = — -det +0 . 3.31 

2(«-«) \u n J' a {u n ) v n J' a (v n )J \ n ) 

In the determinant in ( |3.31| ) we can replace u n and v n by y/u and y/v respectively, and make an 
error which we could estimate using Lemma |3.5| . However, this estimate would not be uniform 
for u — v close to zero. So we will be more careful. We bring in a factor u~ 2 into the first 
column of the determinant in ( |3.31 ) and a factor v~ 2 into the second. Then we subtract the 
second column from the first to obtain 

itfvt , / «~§ J a {u n ) -v~f J a (v n ) v~%J a (v n )\ „fu%v%\ 
D n {u,v) = - -det _ s _„ _ s J+° • 

Z{U-V) y u 2 u n J' a {Un) - V 2V n J' a {Vn) V 2 V n J' a {v n ) J \ U ) 

(3.32) 

From Lemma 3.5 it follows that uniformly for x in bounded subsets of (0, oo), 

— 2 J a (n arccos x n ) — J a {\/x)^ = O ■ 
Then it easily follows that the 1, 1-entry in the determinant in ( |3.32| ) is equal to 
U/u) - v 2 J a U/v) + [ - 



u — V 



n 2 



Similarly, if we use 
d 



dx 



2 n &tccos x n J' a (n aiccos x n ) — x~ ^ y/xJ'^y/x)^ = O (/~2^j > 
we find that the 2, 1-entry is 

U~f y/uj' a {y/u) - y/vj' a (y/v) + O f ~2~ J ■ 



From Lemma 3.5 it also follows that we may replace v n by y/v in the second column at the 
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expense of an error term O(^). Therefore, uniformly for u,v in bounded subsets of (0, oo) 
D n (u,v) 



U 2 V 2 



x det 



2{u - v) 



U 2 



+ o 



a a 



n 



t, , \ uf^S , / « *J a (y/v)-V ™J a (\/v) O (^)\ 

= J a (u, v) + — det a a , . I 

( Q Q \ 

1i 2 t> 2 \ 
—J ' < 3 33 > 

Since z~ a l 2 J a {yfz) is an entire function we get by the mean value theorem that 

u~% J a {Vu) - v~^J a {^/v) 
u — v 

is bounded for u, v in bounded subsets of (0, oo), and similarly, that 
u — V 



is bounded for u,v in bounded subsets of (0, oo). Therefore, we have by ( 3.33| ) 
D n (u, v) = J a (ti, v) + O 



a a 
U2V2 



n I 

uniformly for u, v in bounded subsets of (0, oo), which completes the proof of part (c) of Theorem 
1.1. □ 

3.4 Proof of Corollary 1.2 

To prove Corollary |1.2| we first need two lemmas. 
Lemma 3.6 For u,v S (0, oo) we have as n — > oo, 

The error term is uniform for u, v in bounded subsets of (0, oo). 
Proof. Since 

ufuf / U~* Ja(y/u) - V~% JaWv) V~* J a {y/v) \ 

3 Q {u,v) = — -det a a 

and since v ~ a l 2 J a {y/v) and v ~ a l 2 y/v J' a (\/v) are entire functions of v, it is easy to see, by the 
discussion following fl3.33|), that for every R > 0, there exists a constant c > so that 



S a (u, v) < cu" v™ , for «,t)€ (0, R). (3.35) 
Therefore, ( |3.34| ) is an immediate consequence of Theorem [L.l[ □ 
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In the second lemma, we let D n s and § a s be the integral operators with kernels D n (u,v 
-kzK n (l - 1 - ^i) and J w (u, v) respectively, acting on L 2 (0, s). 

Lemma 3.7 D rhS and are positive trace class operators on L 2 (0,s). 

Proof. Let u n = 1 — ^ and v ra = 1 — ^ • Since 

^ n— 1 

D n (u,v) = — — ^yj w(u n )y/ w(v n ) ^~* j Pj{u n )pj( 



(3.36) 



3=0 



we have that D n)S is a finite rank operator and hence a trace class operator. Now, for every 
/ G L 2 (0, s) we have by ( pOtj ) 



s rs 

J D n (u,v)f(u)f(v)dudv-^^ 



I n ~ l / /•« \2 

^2^2 yj \fw(u n )pj(u n )f(u) duj > 0, (3.37) 



so that D n>s is a positive operator. Letting n — > oo, we have that D n (u,v) — > J a (u, v) for 
every re, t> by Theorem |TTT](c) . By (|H), there is a constant c > independent of re so that for 
re, v G (0, s], 

|D„(«,«)/(it)/(t;)| < cuS«f |/(«)||/(«)|. 

Then by the dominated convergence theorem and ( |3.37 ), 



/ I a (u,v)f(u)f(v)dudv= lim / / D n (u,v)f(u)f(v)dudv>0, 
o Vo n -*°°Jo Jo 

so that J a is positive. Since |J Q (re, v)\ < cu a / 2 v a ^ 2 for re, v G (0, s] we also have 



(3.38) 



tr J QjS | — 


/ S a (u,u)du 


^1 




Jo 


Jo 



which implies that J QjS is a trace class operator. 



□ 



Proof of Corollary 1.2 It is well known that, see for example |3, ffl, 

det (D n (ui, Uk))i< ijk <j dui... duj. 



3=0 



j\ Jo 



For a > 0, we have that D n (u,v) is continuous on [0, s] x [0, s]. Then it follows as in [24| that 
PnO- ~ 1) ^ S e Q ua l to the Fredholm determinant 



(3.39) 



Pn (1-^2,1) = det(I 



For a < 0, the integral kernel D n (u,v) is not continuous, but satisfies an estimate \D n (u,v)\ < 
cu a/2 v a/2 f or u ^ v £ (0, s], for some constant c > 0. This is enough to establish (|3.39j ) also in 
this case. By Theorem 3.4 of p3] we have 



| det (J - £> n , s ) - det (J 



Ill ex P 



where || • is the trace norm in L 2 (0, s). So, in order to prove that det(I — D n s ) converges to 
det(J — I a ,s) as n — ► oo, it is enough to show that Z) niS tends to S a ,s in trace norm. By Theorem 
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2.20 of [24 1 and the positivity of D n s and JJ a s , it then suffices to prove that D n s — > Jq, s weakly, 
and that tr D n ^ s — > trJ QjS as n — > oo. So we have to prove that for f,g € L 2 (0, s), 

lim / / D n (u,v)f(u)g(v)dudv = / / J a {u,v)f{u)g(v)dudv, (3.40) 
n ^°°Jo A) Jo Jo 

and that 

rs rs 

lim / D n (u,u)du= / J Q (ii,ii)du. (3-41) 



Both ( 3.40 ) and ( |3.4lD follow easily from the pointwise convergence D n (u,v) — * J a (u,v), the 



uniform bound \D n (u,v)\ < cu a l 2 v OL l 2 and the dominated convergence theorem. □ 
Acknowledgements 

We thank Ken McLaughlin and Walter Van Assche for useful discussions. 

References 

[1] M. Abramowitz and LA. Stegun, "Handbook of Mathematical Functions," Dover Publica- 
tions, New York, 1968. 

[2] P. Bleher and A. Its, Semiclassical asymptotics of orthogonal polynomials, Riemann-Hilbert 
problem, and universality in the matrix model, Ann. Math. 150 (1999), 185-266. 

[3] P. Deift, "Orthogonal Polynomials and Random Matrices: A Riemann-Hilbert Approach," 
Courant Lecture Notes 3, New York University, 1999. 

[4] P. Deift, T. Kriecherbauer, K. T-R McLaughlin, S. Venakides, and X. Zhou, Uniform 
asymptotics for polynomials orthogonal with respect to varying exponential weights and 
applications to universality questions in random matrix theory. Comm. Pure Appl. Math. 
52 (1999), 1335-1425. 

[5] P. Deift, T. Kriecherbauer, K. T-R McLaughlin, S. Venakides, and X. Zhou, Strong asymp- 
totics of orthogonal polynomials with respect to exponential weights, Comm. Pure Appl. 
Math. 52 (1999), 1491-1552. 

[6] P. Deift, T. Kriecherbauer, K. T-R McLaughlin, S. Venakides, and X. Zhou, A Riemann- 
Hilbert approach to asymptotic questions for orthogonal polynomials, J. Comput. Appl. 
Math. 133 (2001), 47-63. 

[7] P. Deift and X. Zhou, A steepest descent method for oscillatory Riemann-Hilbert problems. 
Asymptotics for the MKdV equation, Ann. Math. 137 (1993), 295-368. 

[8] F. J. Dyson, Correlations between eigenvalues of a random matrix, Comm. Math. Phys. 19 
(1970), 235-250. 

[9] A.S. Fokas, A.R. Its, and A.V. Kitaev, The isomonodromy approach to matrix models in 
2D quantum gravity, Commun. Math. Phys. 147 (1992), 395-430. 

[10] P. J. Forrester, The spectrum edge of random matrix ensembles, Nucl. Phys. B 402 (3) 
(1993), 709-728. 

[11] T. Kriecherbauer and K.T-R McLaughlin, Strong asymptotics of polynomials orthogonal 
with respect to Freud weights, Internat. Math. Res. Notices 1999, no. 6, (1999), 299-324. 



18 



[12] A.B.J. Kuijlaars, K.T-R McLaughlin, W. Van Assche, and M. Vanlessen, The Riemann- 
Hilbert approach to strong asymptotics for orthogonal polynomials, manuscript, 2001. 

[13] H.S. Leff, Class of ensembles in the statistical theory of energy-level spectra, J. Math. Phys. 
5 (1964), 763-768. 

[14] E. Levin and D.S. Lubinsky, " Orthogonal Polynomials for Exponential Weights," CMS 
Books in Mathematics, Vol. 4, Springer- Verlag, New York, 2001. 

[15] D.S. Lubinsky, A survey of general orthogonal polynomials for weights on finite and infinite 
intervals, Acta Appl. Math. 10 (1987), 237-296. 

[16] M.L. Mehta, "Random Matrices," 2nd. ed. Academic Press, Boston, 1991. 

[17] T. Nagao and P.J. Forrester, Asymptotic correlations at the spectrum edge of random 
matrices, Nucl. Phys. B 435(3) (1995), 401-420. 

[18] T. Nagao and K. Slevin, Laguerre ensembles of random matrices: nonuniversal correlation 
functions, J. Math. Phys. 34 (1993), 2317-2330. 

[19] T. Nagao and M. Wadati, Correlation functions of random matrix ensembles related to 
classical orthogonal polynomials, J. Phys. Soc. Japan 60(10) (1991), 3298-3322. 

[20] T. Nagao and M. Wadati, Eigenvalue distribution of random matrices at the spectrum edge 
J. Phys. Soc. Japan 62(11) (1993), 3845-3856. 

[21] P. Nevai, "Orthogonal Polynomials," Memoirs Amer. Math. Soc, Vol. 18, Nr. 213, Amer. 
Math. Soc, Providence R.I., 1979. 

[22] L. Pastur and M. Shcherbina, Universality of the local eigenvalue statistics for a class of 
unitary invariant random matrix ensembles, J. Stat. Phys. 86 (1997), 109-147. 

[23] M. Shiroishi, T. Nagao, and M. Wadati, Level spacing distributions of random matrix 
ensembles, J. Phys. Soc. Japan 62(7) (1993), 2248-2259. 

[24] B. Simon, "Trace Ideals and their Applications," Cambridge University Press, Cambridge, 
1979. 

[25] G. Szego, "Orthogonal Polynomials," Fourth edition, Colloquium Publications, Vol. 23, 
Amer. Math. Soc, Providence R.I., 1975. 

[26] V. Totik, "Weighted Approximation with Varying Weight," Lect. Notes Math. 1569, 
Springer- Verlag, Berlin, 1994. 

[27] C. Tracy and H. Widom, Level spacing distributions and the Airy kernel, Comm. Math. 
Phys. 159 (1994), 151-174. 

[28] C. Tracy and H. Widom, Level spacing distributions and the Bessel kernel, Comm. Math. 
Phys. 161 (1994), 289-309. 



19 



